function warped_points = thinPlateWarp(new_points, control_points, W)

u = thinPlateFunc(sqrt(dist2(new_points, control_points)));
u = [u'; ones(1, size(new_points, 1)); new_points']; 
warped_points = (W * u)';
